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Values  of  broadband  coherence  versus  time  delay  are  presented  for  two  horizontally  separated 
receivers  in  water  overlying  a  rigid  bottom.  These  values  are  given  for  different  source  azimuth 
angles  (from  broadside  toendfire)  and  for  (a)  surface-bottom  (SB)  paths,  (b)  bottom-surface 
(BS)  paths,  and  (c)  a  combination  of  SB  and  bottom-reflected  (B)  paths.  The  acoustic 
wavelength  band  considered  is  the  passband  from  0.2-1. 0  m.  The  water  depth  and  range  are 
both  5  km.  The  source  and  receivers  are  all  at  200  m  depth,  with  the  receivers  separated  by  25 
m.  The  ocean  bottom  is  assumed  to  be  rigid  and  flat.  Scattering  from  the  ocean  surface  is 
modeled  using  the  time-domain,  facet-ensemble  method  [H.  Medwin,  J.  Acoust.  Soc.  Am.  69, 

.  -1060-1064  ( 1981 );  C.  S.  Clay  and  W.  A.  Kinney,  J.  Aroust.  Soc.  Am.  83,  2126-2133  ( 1988) ) 
V  applied  to  a  long-crested  wave  model  of  the  ocean  surface  using  actual  wave-height 
measurements  with  an  rms  value  of  1  m.  A  z-transform  algorithm  (suggested  by  C.  S.  Clayl'is 
used  to  perform  broadband  Altering  and  to  compute  coherence.  Results  demonstrate  how 
signals  that  arrive  endfire  via  SB  paths  possess  greater  coherence  (with  values  near  0.8)  than 
signals  that  arrive  via  BS  paths  (values  near  0).  This  is  because  the  SB  path  involves  scattering 
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off  areas  on  the  ocean  surface  that  are  close  together  and  physically  similar.  /<^ 
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INTRODUCTION 

Two  horizontally  separated  receivers  can  be  used  to  lo¬ 
cate  a  submerged  sound  source  by  correlating  the  received 
signals.  The  signal  arriving  at  each  receiver  does  so  via  multi¬ 
ple  paths  (e  g.,  B,  SB,  and  BS,  see  Fig.  1 ).  The  degree  to 
which  the  signals  at  the  two  receivers  are  correlated  depends 
in  large  part  on  the  physical  features  of  the  boundaries  at  the 
ocean  surface  and  the  ocean  bottom. 

This  article  presents  a  technique  for  estimating  the 
broadband  loss  in  coherence  between  two  received  signals 
due  to  ocean  surface  roughness  alone.  The  technique  utilizes 
the  facet-ensemble  method  for  modeling  rough  surface  scat¬ 
tering.1 2  The  acoustic  wavelength  band  considered  is  from 
0.2-1. 0  m.  In  this  problem,  the  ocean  bottom  is  assumed  to 
be  flat  and  rigid.  Both  the  depth  of  the  water  column  and  the 
range  between  the  source  and  receivers  are  5  km  ( see  Fig.  1 ) . 
The  depths  of  the  source  and  receivers  are  all  200  m  (with 
one  notable  exception  to  be  discussed ) . 

The  problem  of  estimating  the  coherence  loss  at  two 
receivers  due  to  rough  surface  scattering  has  received  pre¬ 
vious  attention  in  the  literature.  Parkins’  formalism3  de¬ 
scribes  the  space  and  time  variation  of  cw  acoustic  signals 
reradiated  from  a  time-varying  ocean  surface  in  terms  of 
local  reflection  from  a  series  of  facets  (i.e.,  the  “broken-mir¬ 
ror”  approach).  Diffraction  effects  at  the  facets  are  neglect¬ 
ed,  and  the  problem  of  broadband  coherence  loss  is  not  ad¬ 
dressed.  McDonald  et  al.*  have  used  the  Fresnel  corrected 
Kirchhotf  approximation  to  compute  the  interfrequency 
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correlation  function  of  scattered  signals.  Their  work  relates 
frequency  shift,  time  lag,  and  spatial  separation  of  the  receiv¬ 
ers.  Again,  the  fully  broadband  correlation  problem  in  the 
time  domain  is  not  considered.  Brill  and  his  co-workers5,6 
have  also  used  a  broken-mirror  approximation  to  calculate 
time,  frequency,  and  angle  spreads  of  acoustic  signals  re¬ 
flecting  from  a  fixed  rough  boundary. 

The  facet -ensemble  method  was  developed  to  overcome 
inaccuracies  in  other  theories  with  regard  to  diffraction,  in 
particular,  and  reflection,  in  general.  The  two  physical  pro¬ 
cesses  are  different  and  cannot  be  treated  in  time-separated 
form  with  the  same  analytical  formalism.  Rough-surface  dif¬ 
fraction  contributions  are  treated  in  the  method  by  using  the 
exact  three-dimensional  solution  for  diffraction  from  a  rigid 
(or  pressure  release )  wedge  or  trough. 1  The  reflected  contri¬ 
butions  from  facets  (wedge  halves)  that  are  specularly  ori¬ 
ented  are  treated  based  on  a  new  interpretation  of  the  Ru- 
binowicz  formulation.2  In  its  present  form,  the 
facet-ensemble  method  calculates  the  three-dimensional 
scattered  field  for  surfaces  possessing  features  for  which 
height  varies  in  only  one  direction  (i.e.,  long-crested  fea¬ 
tures).  The  method  has  been  shown  to  be  accurate  and  ro¬ 
bust  in  comparisons  with  experimental  data  and  is  particu¬ 
larly  well  suited  to  the  problem  being  addressed  in  this 
article.  '■7,K  Currently,  it  is  being  modified  to  handle  three- 
dimensional  surfaces  approximated  using  triangular  facets. 

The  method  can  provide  both  frequency-domain  and 
time-domain-impulse  solutions.  In  this  problem,  a  surface- 
scattered  impulse  response  is  computed  for  each  receiver 
and  then  convolved  (using  a  z-transform  technique)  with  a 
time-domain  broadband  filter.  The  broadband-filtered  sig- 
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nal  at  one  receiver  is  then  correlated  in  time  (again  using  z- 
transforms)  with  the  filtered  signal  at  the  other  receiver.  The 
correlation  is  averaged  over  different  “looks”  at  the  model 
surface.  The  surface  is  deterministically  modeled  in  cross 
section  over  a  spatial  interval  that  is  large  compared  to  the 
ensonified  region  of  interest.  In  the  problem  discussed  here, 
the  surface  cross  section  (a  segment  of  which  is  shown  in 
Fig.  4)  is  taken  from  actual  temporal  waveheight  (wave- 
rider  buoy)  measurements  with  a  root-mean-square  (rms) 
value  of  approximately  1  m.0  The  spatial  waveheight  values 
are  derived  assuming  a  mean  wave  speed  of  7.7  ms'  '.  The 
facet-ensemble  method  uses  a  Monte  Carlo  approach  in  that 
it  averages  over  different  discrete  looks  at  a  deterministic 
representation  of  the  surface. 

Section  I  of  this  article  contains  a  discussion  of  the  spe¬ 
cific  z-transform  techniques  used  to  calculate  the  broadband 
correlation  versus  time  for  two  receivers.  Section  II  provides 
results  for  the  case  of  two  horizontally  separated  receivers  in 
the  presence  of  a  rough  ocean  surface  and  a  flat  rigid  ocean 
bottom.  Results  are  provided  for  individual  cases  involving 
(a)  SB  paths  alone,  (b)  BS  paths  alone,  and  (c)  a  superposi¬ 
tion  of  SB  and  B  paths.  The  results  include  correlation  as  a 
function  of  azimuth  (from  broadside  to  endfire)  and  relative 
depth  of  the  two  receivers.  An  explanation  is  given  as  to  why 
the  correlation  is  better  in  the  case  of  SB  paths  than  it  is  in  the 
case  of  BS  paths.  Section  III  provides  a  set  of  conclusions 
derived  from  the  results. 

I.  THEORY 


/>„(z)  =£a„z"  (n  =  0,1, 2,3,...)  , 
where 

z  =  e  (2) 

where  a>  is  2 n  times  the  sampling  frequency  and  A/  is  the 
sampling  interval,  with  a,„  a ,,  a2 ,  etc.,  being  the  sampled 
signal  amplitudes  at  the  corresponding  time  steps.  A  similar 
series  may  be  written  down  for  the  signal  received  at  B. 

The  correlation  of  the  two  signals  may  formally  be  ex¬ 
pressed  as 

C„„(r)=J  PJt)P*U-r)dt,  (3) 

where  Pa  (/)  and  Ph(t)  are  both  known  over  the  time  win¬ 
dow  (/,->?,)  ofinterest.  In  terms  of  the  z-transform  notation 
of  Eq.  ( 1 ),  the  normalized  correlation  of  the  two  signals  may 
be  defined  as 

Cah(z)  =  Pa(z)P*(z)/N 

=  '£cpzp  (/>  =  o.  ±  1,±2,  ±3,...),  (4) 

P 

where  N  is  the  normalization  factor: 

N  =  (X  l*m|2)  ■  (5) 

'  n  m  / 

The  coefficients  of  Eq.  (4)  are  then  given  by  the  following. 


A  z-transform  technique  is  used  to  perform  the  numeri¬ 
cal  correlation  and  convolution  of  the  time-dependent  sig¬ 
nals  that  will  be  discussed  here.  Since  the  z-transform  meth¬ 
od  is  very  well  established  in  the  literature,10  only  a  brief 
pertinent  outline  of  it  will  be  given. 

Assume  that  an  extended,  time-dependent  signal  P{t), 
arriving  at  receiver  A  in  Fig.  1,  may  be  approximated  by  the 
series 
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CP  =  'Zan+pt>VN  (P>  0), 


cp  =  'Za"bt -p/N  (p<°>- 

n 

The  processing  of  a  signal  (received  at  either  A  or  B) 
through  a  frequency  filter  is  formally  described  by  a  convo¬ 
lution  integral 


?«(')=  r  pa 

Jo 


(r)R(t  —  r)dr , 


where  R  (I)  is  the  impulse  response  of  the  filter,  which  in  z- 
transform  notation  may  be  expressed  as 

R(z)  =  '£rmz”  (m  =  0, 1,2,3,... )  .  (8) 

m 

The  convolution  of  Pa  and  R  is  given  by 

PJz)  =  ^apzf>  (p  =  0,1, 2,3,...),  (9) 

p 

where  the  coefficients  are 


(b) 

TOP  VIEW 


FIG  1  (a)  Schematic  diagram  showing  the  side  view  of  >h*  trice’ 
arrangement  of  the  submerged  source  and  receivers  (A  and  B  at  the  same 
depth  and  separated  by  25  m)  together  with  the  three  types  of  acoustic  ray 
paths  (B,  SB,  BS)  (not  to  scale),  (b)  Top  view  of  the  geometrical  layout 
(with  azimuthal  rotation). 


a„  =  X°"rp  »  («  =  0, 1.2,3 ,...,/>).  (10) 

n 

Computer  algorithms  based  upon  Eqs.  (4)  and  ( 10)  were 
developed  to  perform  the  calculations  discussed  next. 

II.  RESULTS 

The  formalism  discussed  in  the  previous  sections  was 
implemented  to  model  the  correlation  technique  for  the 
problem  shown  schematically  in  Fig.  1 .  An  acoustic  source  is 
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placed  at  a  depth  of  200  m  in  an  isovelocity  ocean  (c=  1500 
ms  1 )  5000  m  deep,  with  a  flat  reflecting  bottom.  Two  hori¬ 
zontally  separated  receivers  spaced  25  m  apart  at  a  depth  of 
200  m  are  placed  at  a  range  of  5000  m  from  the  source. 
Figure  1(a)  shows  that  signals  can  arrive  at  either  receiver 
via  one  of  three  different  types  of  paths:  a  bottom-reflected 
(B)  path,  a  surface-scattered/bottom-refiected  (SB)  path, 
or  a  bottom-reflected/surface-scattered  (BS)  path.  Direct 
paths  and  other  paths  from  the  source  to  the  receivers  that 
involve  multiple  reflections  from  the  ocean  surface  and  bot¬ 
tom  are  not  considered  here.  Figure  1  ( b )  shows  a  top  view  of 
the  geometry.  The  signals  arriving  at  receivers  A  and  B  are 
sampled,  filtered  over  the  passband  of  0.2-  to  1.0-m  wave¬ 
length,  and  correlated  in  time  over  a  range  of  azimuthal  an¬ 
gles  [  0  in  Fig.  1(b)]  from  broadside  (0  =  0  deg )  to  endfire 
(0=90  deg). 

Scattering  from  the  ocean  surface  is  modeled  using  the 
facet-ensemble  method  applied  to  a  long-crested  wave  model 
of  the  surface,  as  described  earlier.  In  the  ocean,  wave  crests 
can  be  typically  30-40  m  long  with  a  crest-to-crest  distance 
of  perhaps  10-20  m.  For  all  of  the  examples  considered  here, 
the  wave  crests  are  assumed  to  be  normal  to  the  direction  of 
propagation.  In  this  regard,  it  is  reasonable  to  assume  that 
the  surface  is  long  crested  across  the  ensonified  regions  for 
the  geometries  considered.  The  normal  azimuthal  orienta¬ 
tion  of  the  crests  was  selected  because  it  provided  the  best 
case  for  correlation  between  the  signals  received,  due  to  the 
fact  that  the  ensonified  regions  are  physically  most  similar. 
Any  other  orientation  of  the  crests  presents  more  physical 
dissimilarity  and  poorer  signal  correlation. 

Figure  2  shows  the  magnitude  of  the  correlation  versus 
time  delay  for  signals  arriving  via  SB  paths  at  the  receivers 
oriented  broadside  and  for  a  perfectly  flat  ocean  surface  ( B 
and  BS  paths  are  discussed  later).  In  this  case,  the  paths 
traversed  are  symmetrically  displaced  on  either  side  of  the 
axis  drawn  from  the  source  to  the  point  midway  between  the 
receivers.  Therefore,  the  two  path  lengths  are  equal,  and  the 
main  correlation  peak  occurs  at  a  correlation  time  of  zero 
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FIG  2.  Correlation  versus  time  delay  for  the  broadside  configuration  (SB 
ray  paths)  and  a  flat  ocean  surface 


TIME  DELAY  (s) 


FIG.  3.  Correlation  curve  for  the  broadside  case  (0  =  Odeg,  SB  ray  paths) 
for  a  rough  ocean  surface  (rms  waveheight  =  1  m) 


(i.e.,  at  the  center  of  the  curve).  Since  the  sea  surface  is 
perfectly  flat,  the  sound  undergoes  specular  reflection,  just 
as  it  does  from  the  ocean  bottom.  Therefore,  the  impulse 
waveform  of  signals  arriving  at  the  receivers  is  identical  to 
that  emitted  by  the  source  (in  this  case  a  narrow  square 
pulse),  but  of  lower  intensity.  After  bandpass  filtering,  the 
main  correlation  peak,  which  has  been  normalized  to  have  a 
maximum  value  of  unity,  has  a  shape  similar  to  a  sine- 
squared  function.  This  can  be  seen  in  Fig.  2.  The  shape  of  this 
correlation  function  is  essentially  determined  by  the  Fourier 
transform  of  the  bandpass  filter.  In  Fig.  3,  the  correlation 
curve  is  shown  for  a  case  with  an  identical  source/receiver 
configuration  to  that  used  to  obtain  Fig.  2.  Now,  however, 
the  sea  surface  is  no  longer  flat,  and  the  surface  waves  are 
modeled  in  cross  section  and  with  an  rms  waveheight  value 
of  1  m.  Figure  4  shows  a  sample  of  the  surface  cross  section 
which  was  derived  using  the  wave-rider  buoy  data  men¬ 
tioned  earlier.9  The  wave  crests  are  normal  to  the  direction 
of  propagation  and.  as  before,  only  the  SB  path  is  considered. 
The  path  lengths  from  the  source  to  the  two  receivers  are 
again  identical,  so  that  a  main  peak  is  observed  at  a  correla¬ 
tion  time  of  zero.  However,  since  the  sea  surface  is  not 
smooth,  the  scattered  waves  have  an  impulse  waveform  re- 


0  50  100  150  200  250  300  350  400  450 


HORIZONTAL  DISTANCE  (m) 

FIG.  4  Cross  section  of  portion  of  ocean  surface  obtained  from  wave-rider 
buoy  data.  This  surface  ( rms  height  -  1  m )  was  used  to  derive  the  correla¬ 
tion  plots. 


1032 


J.  Acoust  Soc  Am  .  Vol  84.  No  3,  September  1988 


C  Feumade  and  W  A  Kinney:  Broadband  temporal  coherence 


1032 


suiting  from  the  aggregation  of  waves  scattered  from  various 
facets  on  the  sea  surface.  Since  the  signals  are  not  simply 
specularly  reflected  from  a  flat  sea  surface,  the  impulse 
waveforms  arriving  at  the  receivers  are  radically  different 
from  the  narrow  spike  emitted  by  the  source.  Therefore,  the 
main  correlation  peak  (which  is  again  symmetric  around 
t  =  0  and  has  a  maximum  value  of  unity)  no  longer  has  the 
simple  shape  seen  in  Fig.  2.  Instead,  there  is  significant 
broadening  at  the  base  of  the  peak  and  considerable  random 
variation  stretching  out  both  before  and  after  the  main  peak. 
This  is  due  to  the  cross  correlation  of  delayed  signals  scat¬ 
tered  from  different  portions  of  the  ocean  surface. 

Figure  5  demonstrates  the  change  in  the  correlation  co¬ 
efficient  as  the  azimuthal  angle  of  the  source  is  increased 
from  0  deg  (broadside)  to  90  deg  (endfire).  Again,  the 
wavecrests  are  normal  to  the  propagation  direction  and  only 
SB  paths  are  considered.  Two  new  features  appear  in  this 
figure.  First,  although  the  curve  is  again  normalized  to  unity 
as  defined  by  the  peak  correlation  coefficient  calculated  in 
the  broadside  case  (see  Fig.  3),  the  maximum  correlation 
coefficient  value  has  decreased  from  unity  to  about  0.55  due 
to  the  azimuthal  rotation.  The  reason  for  this  is  that  the 
signals  reaching  the  two  receivers  are  no  longer  identical,  but 
differ  because  they  have  been  scattered  from  different  se¬ 
quences  of  surface  waves.  Not  only  is  the  peak  value  re¬ 
duced,  but  the  symmetry  of  the  shape  of  the  correlation 
curve  about  the  peak  (observed  in  Figs.  2  and  3 )  is  also  lost. 
Second,  there  is  a  shift  in  the  position  of  the  correlation  peak 
away  from  /  =  0  The  change  in  the  azimuthal  angle  causes  a 
corresponding  change  in  the  acoustic  path  lengths  from  the 
source  to  the  receivers.  One  receiver  moves  towards  the 
source  (receiver  B  in  Fig.  1),  while  the  other  (receiver  A) 
moves  away.  The  corresponding  difference  in  the  two  travel 
times  gives  the  shift  in  the  peak  position.  In  this  case,  its 
position  corresponds  to  a  time  delay  of 0.0075  s  between  the 
arrivals  of  the  signals  at  receivers  A  and  B.It  should  be  noted 
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FIO  5  Correlation  curve  for  the  endfire  case  (0  =  90  deg.  SB  paths)  and 
for  a  rough  ocean  surface  ( rms  waveheight  I  m ) 


FIG.  6.  Correlation  curve  for  the  broadside  case  (0  =  0  deg)  for  a  rough 
ocean  surface  ( rms  waveheight  =  1  m ) ,  using  BS  paths  only. 


that  even  for  this  extreme  endfire  case  the  correlation  is  still 
good. 

Figures  6  and  7  display  a  set  of  correlation  curves  simi¬ 
lar  to  those  of  Figs.  3  and  5  except  that  now  only  bottom 
surface  ( BS)  paths  from  the  source  to  the  receivers  are  con¬ 
sidered,  rather  than  SB  paths.  In  Fig.  6,  the  broadside  case  is 
presented.  Here,  again,  there  exists  perfect  correlation:  a 
maximum  value  of  unity  at  a  time  delay  of  0  s.  This  curve  is 
similar  to  that  of  Fig.  3  for  the  SB  case.  However,  marked 
differences  between  the  SB  and  BS  cases  begin  to  appear 
when  the  azimuthal  angle  is  increased  away  from  broadside. 
This  is  illustrated  in  Fig.  7,  which  shows  the  correlation 
curve  for  an  azimuthal  angle  of  only  3  deg  for  the  BS  paths. 
Figure  7  should  be  compared  with  Fig.  5,  which  shows  the 
corresponding  curve  for  the  SB  case  for  endfire.  Whereas,  in 


FIG  7  Correlation  curve  for  an  azimuthal  rotation  of  0  -  3  deg  from 
broadside,  using  BS  paths  only 
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the  SB  case,  the  peak  correlation  coefficient  drops  from  uni¬ 
ty  to  near  0.55  for  an  azimuthal  angle  of  90  deg,  in  the  BS 
case,  it  is  approximately  0. 1  for  an  angle  of  only  3  deg. 

The  rapid  fall  in  the  amplitude  of  the  correlation  peak  in 
the  BS  case  for  increasing  azimuth  angle  is  best  explained  by 
examining  the  different  paths  of  propagation.  In  the  SB  case, 
Fig.  8  shows  in  schematic  form  the  paths  taken  by  the  acous¬ 
tic  rays  as  they  travel  from  the  source  to  the  two  receivers  in 
the  broadside  configuration.  The  upper  diagram  shows  a 
side  view  of  the  environmental  geometry.  Note  that  the 
shortest  SB  path  between  the  source  and  either  of  the  two 
receivers  is  taken  by  the  ray  that  is  specularly  scattered  from 
both  the  sea-surface  and  the  sea-bottom  interfaces.  This  is 
denoted  the  “median”  path  in  the  figure.  Rays  traveling 
along  this  path  take  the  least  time  to  pass  from  source  to 
receivers  via  the  surface  and  bottom.  Rays  that  deviate  from 
the  median  path  take  longer  to  pass  from  the  source  to  the 
receivers.  A  time  interval  can  be  specified  such  that  the  earli¬ 
est  median  path  signals  arrive  at  the  inception  of  the  time 
interval.  The  end  of  the  time  interval  can  then  be  used  to 
define  two  “extreme”  paths  (see,  again,  Fig.  8).  The  points 
on  the  sea  surface  at  which  these  two  extreme  paths  reflect 
provide  the  range  limits  of  two  regions  of  interaction 
between  the  rays  and  the  surface.  Sound  scattered  from  these 
regions  may  reach  the  detectors  within  the  specified  time 
period.  The  amount  of  energy  arriving  from  outside  these 
regions  is  small  and  can  be  neglected.  The  second  of  the 
diagrams  in  Fig.  8  shows  the  top  view  of  the  environment. 
Physically,  scattering  from  the  sea  surface  will  occur  from 
areas  surrounding  the  points  where  the  median  paths  inter¬ 
sect  the  surface.  This  is  represented  schematically  in  Fig.  8 
by  the  shaded  elliptical  areas.  The  signal  processing  proce¬ 
dure  used  here  consists  of  performing  phase  and  amplitude 
comparisons  of  signals  received  at  A  and  B.  Since  the  two 
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FIG.  8.  Ray  paths  and  regions  of  acoustic  interaction  with  the  ocean  sur¬ 
face,  for  the  SB  case  in  the  broadside  configuration  ( 6  =  Odeg):  (a)  side 
view,  (b)  top  view 


signals  are  scattered  from  different  regions  of  the  surface, 
this  involves  the  assumption  that  there  is  some  degree  of 
phase  correlation  or  coherence  for  acoustic  waves  scattering 
from  different  but  closely  adjacent  regions  of  the  sea  surface. 
The  degree  of  coherence  is  also  dependent  upon  the  wave- 
heights  and  randomness  of  the  ocean  surface.  In  general, 
however,  the  greater  the  separation  between  the  two  regions 
on  the  surface,  the  less  correlation  there  will  be. 

In  the  BS  case,  an  envelope  containing  the  rays  that 
carry  most  of  the  acoustic  energy  between  source  and  receiv¬ 
ers  may  be  constructed  in  a  similar  manner  to  the  SB  case. 
Again,  a  specified  time  interval  defines  the  median  and  ex¬ 
treme  paths  that  form  the  ray  envelope.  However,  in  the  BS 
case,  there  is  an  important  difference.  In  the  BS  case,  the  rays 
interact  with  the  surface  after  they  are  reflected  from  the 
bottom,  such  that  the  regions  of  ocean  surface  interaction 
are  situated  closer  to  the  receivers  than  to  the  source.  (This 
may  be  visualized  by  mentally  interchanging  the  receiver 
position  with  the  source  position  in  Fig.  8(a)  and  reversing 
the  direction  of  the  arrowheads.  ]  The  overall  size  of  the  re¬ 
gions  of  interaction  is  the  same  as  in  the  SB  case  (determined 
by  the  chosen  time  window,  as  described  above),  but  they 
are  much  more  widely  separated  laterally  from  each  other 
due  to  their  closer  proximity  to  the  receivers.  The  difference 
in  the  correlation  effects  between  the  SB  and  BS  cases  is  not 
noticeable  at  broadside  because,  in  both  cases,  the  two  sig¬ 
nals  travel  to  their  respective  receivers  along  effectively  iden¬ 
tical  paths  due  to  the  long-crested  model  of  the  ocean  surface 
that  has  been  adopted.  This  means  that  perfect  correlation  is 
observed  at  zero  time  delay.  However,  significant  differences 
are  observed  as  the  azimuthal  angle  deviates  from  broadside. 
In  Fig.  9,  the  effect  of  an  azimuthal  rotation  of  the  receivers 
on  the  positions  of  the  regions  of  interaction  is  demonstrated 
for  the  SB  and  BS  cases.  It  is  clear  from  the  diagram  that,  due 
to  the  closer  proximity  of  these  regions  to  receivers  A  and  B 
in  the  BS  case,  they  are  laterally  displaced  with  respect  to 
each  other  much  more  in  the  BS  than  in  the  SB  case.  This 
means  that  sound  scattered  from  these  regions  onto  their 
respective  receivers  will  have  been  scattered  from  surface 


FIG.  9.  Top  views  of  the  regions  of  acoustic  interaction  for  the  SB  and  BS 
cases  after  azimuthal  rotation 
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wave  sequences,  within  the  regions  of  interaction,  that  are 
less  correlated  physically  with  each  other  in  the  BS  case  than 
in  the  SB  case.  In  other  words,  the  two  regions  of  interaction 
are  more  dissimilar  in  the  BS  case,  and  this  gives  rise  to  the 
rapid  fall  in  correlation  with  increasing  azimuth  angle. 

An  effect  similar  to  that  induced  by  azimuthal  receiver 
rotation  is  observed  when  the  receivers  are  held  in  the  broad¬ 
side  configuration,  but  the  depth  of  one  of  them  is  slightly 
increased  with  respect  to  the  other.  In  Fig.  10,  the  correla¬ 
tion  curve  for  the  SB  case  is  shown  where  receiver  A  is  held 
1.5  m  deeper  than  receiver  B.  Comparison  of  this  figure  with 
Fig.  3  shows  that  the  effect  is  to  decrease  the  amplitude  of  the 
correlation  peak  to  0.9,  and  to  move  its  position  slightly 
away  from  zero  time-delay.  Figure  1 1  shows  the  correspond¬ 
ing  curve  for  the  BS  case  where  a  dramatically  decreased 
correlation  peak  amplitude  is  seen.  The  explanation  is  much 
the  same  as  for  rotation  in  azimuth.  Dropping  one  receiver 
by  1.5  m  introduces  a  much  greater  difference  in  the  surface 
wave-profile  encountered  by  the  acoustic  rays  in  the  BS  case 
than  in  the  SB  case.  The  signals  reaching  the  two  receivers 
will  therefore  be  less  correlated  in  the  BS  case. 

Another  interesting  difference  between  Figs.  10  and  1 1 
should  also  be  noted.  In  Fig.  10,  the  correlation  peak  has 
been  moved  to  a  positive  time  delay  from  zero;  whereas,  in 
Fig.  1 1,  it  has  been  moved  to  a  negative  time  delay.  This  is 
because,  when  receiver  A  is  held  1.5  m  deeper  than  receiver 
B,  receiver  A  is  reached  first  by  the  rays  in  the  SB  case, 
whereas  it  is  reached  second  in  the  BS  case. 

All  of  the  correlation  curves  shown  have  been  formed  by 
the  comparison  of  two  signals  that  are  of  the  same  type  (ei¬ 
ther  both  SB  or  both  BS).  Figure  12  shows  the  correlation 
curve  for  signals  arriving  at  receivers  A  and  B  which  are 
comprised  of  rays  that  have  traveled  along  both  SB  and  B 
paths.  It  must  be  remembered  that  the  time  response  of  the 
bottom-reflected  rays  will  be  a  short  square-wave  spike  emit¬ 
ted  by  the  source,  but  decreased  in  amplitude  due  to  spheri- 


FICi.  10  Correlation  for  the  broadside  case  (SB  paths)  for  a  rough  <x*can 
surface  (rms  waveheight  -  1m),  and  with  one  receiver  held  15  m  deeper 
than  the  other 
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FIG  1 1  Correlation  curve  for  the  broadside  case  ( 6  =  Odeg,  BS  paths),  for 
a  rough  ocean  surface  ( rms  waveheight  =  1  m ) ,  and  with  one  receiver  held 
1.5  m  deeper  than  the  other. 


cal  spreading.  This  is  because  reflection  from  the  bottom 
(which  is  flat  and  rigid)  is  taken  here  to  be  simply  specular 
and  introduces  no  incoherent  scattering.  In  Fig.  12,  the  cor¬ 
relation  curve  is  shown  for  the  two  receivers  at  the  same 
depth  in  the  broadside  configuration  and  for  a  perfectly  flat 
sea  surface.  At  the  position  of  zero  time  delay,  a  central  cor¬ 
relation  peak  appears  that  is  formed  by  the  addition  of  the 
correlation  peak  for  the  two  SB  signals  together  with  the 
corresponding  peak  for  the  two  B  signals  together.  In  addi¬ 
tion,  two  other  peaks  appear  at  time  delays  of  +  0.24  s. 
These  peaks  are  caused  by  the  correlation  of  the  SB  signal  at 


TIME  DELAY  (s) 


FIG  12.  Correlation  for  the  broadside  case  ( 0  -  0  deg)  with  a  flat  ocean 
surface  Thecioss  correlation  is  shown  between  both  SB  and  B  path  signals 
T  he  high  central  peak  is  due  to  SB-SB  and  B- B  cross  correlations  The  two 
smaller  side  lobes  are  due  to  SB  B  and  B-SB  c-nv,  correlations 
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receiver  A  with  the  B  signal  at  receiver  B,  and  vice  versa.  The 
main  peak  has  a  normalized  value  of  1  and  the  side  peaks 
have  values  of  0.5  each.  These  are  the  values  one  would  ex¬ 
pect  for  a  perfectly  flat  ocean  surface  and  ocean  bottom.  For 
the  same  configuration  of  source  and  receiver  positions,  but 
with  a  rough  sea  surface  with  waves  of  rms  height  1  m,  the 
level  of  the  side  peaks  drops  from  0.5  to  0.21.  There  is  no 
appreciable  fall  in  the  level  of  the  central  peak  (cf.  Figs.  2 
and  3).  For  the  endfire  configuration  (with  rough  surface), 
which  is  the  worst  geometrically,  the  amplitude  of  the  cen¬ 
tral  peak  falls  to  0.63  but  is  still  clearly  distinguishable  from 
the  sidelobes,  which  have  an  average  value  ==0.3.  Also,  even 
though  the  central  peak  is  now  formed  by  the  addition  of  the 
correlation  peak  for  the  two  SB  signals  together  with  the 
peak  for  the  two  B  signals  together,  these  two  peaks  do  not 
separate  in  the  endfire  configuration.  Therefore,  informa¬ 
tion  contained  in  the  main  peak  should  be  usable. 

III.  CONCLUSIONS 

A  z-transform  technique  has  been  presented  and  used  in 
conjunction  with  the  time-domain  facet-ensemble  method  to 
compute  coherence  versus  time  delay  for  two  receivers.  It 
showed  that,  for  a  rough  ocean  surface  and  a  flat  rigid  bot¬ 
tom,  signals  arriving  via  SB  paths  possess  greater  coherence 
than  signals  arriving  via  BS  paths.  The  dissimilarity  between 
the  two  scattering  regions  on  the  ocean  surface  for  the  two 
receivers  is  greater  in  the  case  of  the  BS  paths  than  the  SB 
paths.  This  implies  that,  for  the  purpose  of  correlating  the 
signals  arriving  at  the  two  receivers  in  the  presence  of  ocean 
surface  roughness  (with  a  flat  rigid  bottom),  it  is  more  ad¬ 
vantageous  to  use  bottom  arriving  energy  than  surface  arriv¬ 
ing  energy.  This  remains  true  when  signals  arriving  via  SB 
and  B  paths  are  superposed. 
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